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ABSTRACT 


An experimental study of the turbulent mixing of subsonic axisym- 
metric gas streams was conducted. Hydrogen-air and air-air mixing 
systems were studied, and the velocity ratio (jet velocity/outer stream 
‘velocity) was varied from 2.4 to 6.3. Special emphasis is placed on 
(1) the centerline decay and radial profile shapes of composition, velocity, 
and total enthalpy and (2) the relationships between the turbulent trans- 
port of mass, momentum, and energy. The major conclusions drawn for 
this particular set of*conditions are (1) for the hydrogen-air system the 
centerline decay decreases with increasing velocity ratio, (2) profile 
similarity of composition, velocity, and total enthalpy is a valid assump- 
tion, (3) there is a definite relationship between the transport of momen- 
tum and energy which is not compatible with a constant Prandtl number, 
and (4) unity Lewis number is a valid assumption. The experimental 
data obtained are tabulated for the benefit of other investigators, 
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NOMENCLATURE 
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Hydrogen mass fraction 
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Inner nozzle diameter 
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Total temperature in the inner plenum 

Velocity in the axial direction 


Maximum velocity in Prandtl eddy viscosity model 
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Umin Minimum velocity in Prandtl eddy viscosity model 

Vv Velocity in the radial direction 

W Hydrogen mass flow rate from mass balance 

W; Hydrogen mass flow rate from metering orifice measurement 
x Axial distance 

Xo Potential core length 

Xoc Potential core length from composition 

Xou Potential core length from velocity 

Y Ratio of specific heats 

A Combined boundary-layer thickness 


Eddy viscosity 


r Mass flux ratio, pjU;/PQUo 
p Density 
¢ H = Hy 
H Hj - Ho 
u-u 
re) 
%u Uj - Uo 
u-u 
Vy : 
Ue aa Ug 
W Molecular weight 
SUBSCRIPTS 
c Centerline 
j Jet 
0) Outer stream 
s Static 
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SECTION I 
INTRODUCTION 


The turbulent mixing between high-speed coaxial streams is of 
primary importance in many engineering devices, such as jet pumps, 
ejectors, and ramjet combustors. The flow in these devices must at 
best-be analyzed by semiempirical techniques because there is no funda- 
mental and complete theory for turbulent flows. Numerous investiga- 
tions of the turbulent mixing process have been conducted, but they have 
failed to produce a generalized theory. The main reason for the lack of 
success in solving this problem is that the turbulent transport proper - 
ties are a function of the fluid dynamics of the flow system. 


The case of a single jet mixing with a quiescent medium has been 
given considerable attention by a number of investigators. Most of 
these studies were limited to incompressible flow, and the results are 
well covered by Schlichting (Ref. 1) and Pai(Ref. 2). One of the first 
investigations of the case of two coaxial streams mixing together was 
conducted by Forstall and Shapiro (Ref. 3). An inner stream composed 
of air with ten percent by volume of helium as a tracer was mixed with 
an outer stream of air, The gas velocities were in the low subsonic 
range, and the temperatures were maintained nearly equal. The main 
conclusions were that (1) momentum is transported less rapidly than 
mass and (2) the normalized velocity and composition profiles exhibit 
shape similarity; that is, the nondimensional radial profiles are in- 
variant with axial distance. 


Alpinieri (Ref. 4) obtained experimental data on the turbulent mixing 
between carbon dioxide and hydrogen central jets exhausting into a 
moving concentric stream of air. The flow velocities were in the low to 
high subsonic range, and the temperatures of the streams were approxi- 
mately equal. Radial and axial distributions of composition and velocity 
were presented. Alpinieri verified Forstall and Shapiro's conclusion 
that mass is transported more readily than momentum. Also, Alpinieri 
concluded that the mixing rate does not approach zero when either the 
velocities or the mass fluxes of the stream are equal. Conclusions con- 
trary to this were suggested by turbulent eddy viscosity models proposed 
by other investigators such as Prandtl as presented by Schlichting (Ref. 1) 
and Ferri (Ref. 5). The eddy viscosity models in question are as follows: 


1. Eddy viscosity model proposed by Prandtl 


€= kb (umax - Umin) 
2, Eddy viscosity model proposed by Ferri 


ea= © (poo = Pole) 
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Zakkay and others (Ref. 6) conducted an experimental investigation 
to determine the turbulent transport coefficients for hydrogen-, helium-, 
and argon-air mixing systems. Central jets of hydrogen, helium, and 
argon at subsonic velocities were injected into an outer stream of air 
maintained at a constant Mach number of 1.6. The ratio of the inner jet 
velocity to the outer stream velocity for the hydrogen-air mixing sys- 
tem was varied from 0.768 to 2.42. In addition to turbulent transport 
coefficients, centerline decay of velocity and composition was presented. 
A comparison between Zakkay's results and results from the current in- 
vestigation is given in Section 3. 3. 


1.1 OBJECTIVE AND JUSTIFICATION 


The objective of this experimental investigation was to document 
the turbulent mixing of subsonic axisymmetric hydrogen and air streams 
at velocity ratios which have not been previously reported. This is part 
of a more extensive investigation of turbulent mixing systems. The ob- 
jective of the overall investigation is to develop empirical data on free 
turbulent flows over a broad range of conditions, so that more percep- 
tive semiempirical theories can be developed. Special emphasis is 
placed on (1) determining the effect of velocity ratio and density gradi- 
ents on the centerline decay of composition, velocity, and total enthalpy, 
and on the composition, velocity, and total enthalpy profile shapes, and 
(2) presenting nondimensional composition, velocity, and total enthalpy 
relationships which are indicative of the turbulent Prandtl, Lewis, and 
Schmidt numbers. 


Integral techniques for solving turbulent mixing problems (Refs. 7 
and 8) employ the assumption that composition, velocity, and total en- 
thalpy profiles exhibit shape similarity. Forstall and Shapiro (Ref. 3). 
listed the following relationships as being representative of their velocity 
profile shapes: 


Cosine curve 


u=% 4 (1 7R 
Ug — Uy 2 , ve 2 Rau 


) 


Three-halves power curve 


2 
=. [1 — 0.293 (R/Rmu)*/ 4] 


Ue ~ Ug 


and error curve 
(R/Ray)? 


ae Le] 
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If‘ profile shape similarity is a valid assumption, a representative pro- 
file ‘curve, the centerline and outer stream velocities, and a character - 
istic width are sufficient to determine the radial velocity distribution at 
any axial location. 


“In most proposed methods for solving turbulent mixing problems, 
the ‘turbulent Prandtl and Lewis numbers are assumed to be unity to 
simplify the mathematical procedure. If the Prandtl and Lewis num- 
bers are unity, the Schmidt number must be unity by definition. If the 
Prandtl and Lewis numbers are unity and the laminar transport coeffi- 
cients are replaced by the corresponding turbulent values, the boundary- 
layer equations for axisymmetric flow may be written as follows: 


Momentum equation, © 


a gu Pi: ae ae Gu) _ OP 
ead pu ox * PY oR 1 Pilpe rR :) x (1) 
Energy equation, 
oH . a 1 @ oH 
pe He py Mad B(en gH) (2) 
Conservation of elemental species, 
gC GA, ae gc 
. pe + pv = 7 Bloen Se) (3) 
and 
Global continuity equation, 
0 (pu) 1 @ (pyR) 
’ pacat Sees or IO 
Ox . R OR (4) 


For the case of constant pressure mixing of an initially uniform infinite 
‘stream, the 8P/@x term in Eq. (1) is zero, and Eqs. (1), (2), and (3) 
are identical in form. If the reference conditions are constant and the 
boundary conditions are similar, a linear relation may be obtained be- - 
tween the variables u, H, and C: 


Uj — Uo . H;— Ho C;- Cy (5) 


The experimental relationships between the terms of Eq. (5) are pre- 
sented in this investigation as being indicative of the Prandtl and Lewis 
number variation. 
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1.2 APPROACH 


In general, turbulent mixing is influenced by the following factors: 
1. Velocity ratio between the streams, 
2. Density gradients in the mixing region, 
3. Axial and radial static pressure gradients, and 
4. Initial boundary-layer and free-stream turbulence level. 


The approach which has been taken in this investigation is to reduce the 
effect of initial boundary-layer and static pressure gradients and to con- 
centrate on the effects of velocity ratio and the density field. The 
initial boundary-layer effect was reduced by designing the nozzle to 
minimize the boundary-layer thickness at the entrance to the test 
section. The initial boundary-layer results are discussed in Section 3.6. 
The static pressure variation was minimized by exhausting to atmos- 
phere as a free jet and maintaining the flows subsonic so that shock 
waves were not formed. However, the static pressure in the mixing 
region was measured, and the gradients near the nozzle exit were 
larger than anticipated. Static pressure data are presented in 

Section 3. 7. 


Hydrogen-air and air-air turbulent mixing systems were investi- 
gated. The hydrogen-air system was chosen to provide a system with 
very large density gradients. Also, the high speed of sound of hydrogen 
makes it possible to attain high velocities while the Mach numbers re- 
main subsonic. Furthermore, hydrogen-air mixtures are frequently 
used in combustion processes, and the results should be useful when 
chemically reactive systems are investigated. The air-air system is 
studied because it provides a system with small density gradients to 
use for comparison with the high density gradient system. 


The approach of this investigation required that a large amount of 
experimental data be obtained. The experimental measurements are 
time-mean average values because techniques for measuring the 
fluctuating quantities are not well developed, especially in the high 
velocity region. 


SECTION II 
APPARATUS AND PROCEDURE 
2.1 FREE-JET MIXING TEST CELL 


A schematic diagram of the free-jet mixing test cell is shown in 
Fig. 1 (AppendixI). Air, which may be heated to 1500°R by an 
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indirect-fired heater, flows around the inner plenum and nozzle. It then 
passes through a 3.5-in. -diam subsonic nozzle to form an annulus 
around the subsonic flow from the inner nozzle. The inner nozzle diam- 
eter is 0.5 in., and the thickness of the trailing edge is 0.005 in. The 
inner nozzle and outer nozzle configuration was designed to produce as 
small an initial boundary layer at the entrance to the test section as 
practical. The inner and outer nozzles were aligned to give flow with 
centerlines which are parallel within less than 0.5 deg. The alignment 
was checked by means of total pressure measurements in the down- 
stream flow field. The test section is open to the atmosphere, but the 
gases from the nozzle are removed by a downstream scoop attached to 
the RTF exhaust system. 


2.2 INSTRUMENTATION 


A Systems Engineering Laboratory (SEL) 600 data acquisition sys - 
tem was used to record all of the data in millivolts on magnetic tape. 
The temperatures were measured with copper-constantan and iron- 
constantan thermocouples, and the pressures were measured with 
strain-gage-type transducers. The gas composition was measured 
with a fluid oscillator, which was developed by the RTF Research 
Branch (Ref. 9). The probe positions were indicated by wire-wound 
potentiometers mechanically connected to the probe, and the inner 
stream flow rate was measured with a calibrated choked orifice. Esti- 
mates of the accuracy of the measured parameters are presented in 
Appendix III. 


A dual-probe arrangement was used to measure the total pressure, 
total temperature, gas compoSition, and static pressure at various 
stations throughout the flow field. A photograph and a sketch of the 
dual-probe arrangement are shown in Fig. 2, and a schematic diagram 
of the probe-related components is shown in Fig. 3. The probe used 
to measure total pressure, total temperature, and gas composition is 
operated in two modes. ‘Total temperature and gas composition are 
recorded on one mode when the probe is aspirated to a vacuum source, 
and total pressure is recorded on the other mode when there is no flow 
through the probe. Static pressure is recorded during both modes of 
operation. The static pressure measurements as recorded are dis- 
placed 0.5 in. to one side of the total pressure probe measurements. 
They are shifted in the data reduction program to align with the total 
pressure probe measurements. 


The dual-probe arrangement is actuated in the flow field by a three- 
position probe actuatoz. The probe location and the initial test condi- 
tions were monitored on an X-Y-Y plotter and strip-chart recorders. 
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2.3 TESTING PROCEDURE 


In-place calibrations were made on all of the pressure transducers 
and potentiometers before each test. The temperature channels were 
calibrated by applying two different millivolt levels. Next, the test 
conditions were established by setting the total pressure in the inner 
and outer plenums and the total temperature of the outer airstream. 
The conditions were allowed to stabilize, and then the probe was set 
on the vertical centerline of the flow field by the following procedure: 
The probe was located on the approximate horizontal centerline of the 
jet by using the peak of the horizontal total pressure profile. Then the 
probe was moved axially to the decaying region of the jet and actuated 
vertically. The vertical centerline was taken to be the location of 
maximum pressure if the jet total pressure was greater than the outer 
stream total pressure or the minimum pressure if the jet total pres- 
sure was below the outer stream pressure. 


The data were recorded by a data acquisition system which is oper- 
ated in the following manner: The probe is traversed radially at a fixed 
axial location until the desired radial location is reached. A single 
switch starts a sequence which stops the probe and begins to record 
each data item 300 times per second. The probe is put into the total 
temperature and gas composition mode of operation, then into the total 
pressure mode. The time in each mode of operation is indicated by a 
function signal switch recorded on magnetic tape. The magnetic tape 
drive is stopped automatically when all data at the radial location have 
been recorded. The probe drive is engaged automatically to traverse 
the probe to the next radial data point location selected by the test con- 
ductor. The total time required to obtain a data point in a radial sur- 
vey is approximately 5 sec, After data points have been recorded to 
give complete profiles at a given axial location, the radial traverse of 
the probe is stopped, and the probe is moved to another axial location 
where the next radial profile is obtained. 


2.4 DATA REDUCTION PROCEDURE 


The data are reduced in three steps with the aid of a digital com- 
puter. First, the magnetic tape from the SEL 600 tape system is pro- 
cessed through a data reduction program which uses the calibrations to 
convert the millivolt signals to engineering units of pressure, tempera- 
ture, and probe position. The data are averaged over 0, 167-sec in- 
tervals (50 scans) to obtain the mean values. The data recorded during 
stabilization of the parameters were discarded. A printout and a tape 
are made of the resulting engineering units data. 
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Second, the engineering units tape is reduced by computer to give 
gas mixture properties at each probe position. The properties obtained 
are composition, density, velocity, and total enthalpy at the initial 
stream conditions and at local points in the flow field. The methods of 
calculating the specific properties are given in Appendix IV. These 
basic properties are recorded on the tape to be used as the inputs for 
the final reduction program. The method by which the data are further 
reduced may change as the knowledge of turbulent mixing increases. 
For this reason, the basic properties tape is stored so that other 
methods of data reduction may be applied to it in the future. 


Finally, the basic properties tape is reduced to provide profiles of 
nondimensionalized velocity, composition, and total enthalpy as well as 
relationships between the transport of mass, momentum, and energy. 
Also, a hydrogen mass balance is made at each axial station recorded 
by integrating the hydrogen mass flux (puC) over the radial distance to 
a control volume outside the mixing zone. The mass balance is com- 
pared with the hydrogen flow measured with a calibrated choked orifice 
to give a consistency check. The results of the consistency check are 
given in Appendix V. 


SECTION III 
DISCUSSION OF RESULTS 


3.1 FLOW FIELD DESCRIPTION 


‘The jet gas and the outer stream gas mix in the inner mixing zone 
as shown schematically in Fig. 4. For all cases investigated, the jet 
gas used was either hydrogen or air at ambient temperature. Only 
regimes I and II are considered in this investigation. In regime III, the 
conditions in the outer stream are a function of the outer stream mixing 
with the surrounding quiescent air. This provides a problem with dif- 
ferent boundary conditions from the problem under consideration. The 
inner potential core length varies with the test conditions. The inner 
potential core is defined as the region in which the composition and 
velocity are constant and equal to the inner nozzle exit conditions. 


3.2 EXPERIMENTAL DATA 


‘Experimental data were obtained for nine different test conditions 
The velocities, total pressures, and total temperatures for these tests 
are tabulated in Table I (Appendix II). These data are divided into three 
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test series. Series I and II are hydrogen-air mixing tests. For the 
first five test conditions (Series I), the total pressure of the two streams 
was set to give a range of velocity ratios from 2.4 to 6.3. In Series I, 
the next three test conditions, the outer stream temperature was raised 
from 650 to 1050°R. The total pressures were varied to give velocity 
ratios from 2.5 to 4.6. Series III, the last test condition, was an air- 
air mixing test at a velocity ratio of 2.4. 


In the near field region* the gradients of composition, total pres- 
sure, static pressure, and total temperature were very large. Insuffi- 
cient data were obtained to define the profiles adequately because the 
increment traveled between radial points was too large. For this 
reason, much of the data in this region is omitted from the results 
presented. The main difficulty in using the near field data is that the 
radial centerline of the flow cannot be accurately determined. The 
experimental centerline is determined by fitting the center portion of 
the composition distribution with an exponential curve and using the 
center of the exponential curve as the centerline of the flow field. When 
there are insufficient data to define the curve, this procedure gives 
centerlines which are obviously in error. 


The experimental data are presented in tabular form in Appendix VI. 
A discussion of the probable accuracy of the experimental measure- 
ments is presented in Appendix III. 


3.3 CENTERLINE DECAY 


The centerline decay of composition and velocity for all of the 
hydrogen-air mixing tests is shown in Figs. 5 and 6, respectively. The 
composition and the velocity are nondimensionalized in the same manner. 
However, for the hydrogen-air system, Cj = 1 and Co, = 0; therefore, 


These curves indicate that the centerline decay decreases with in- 
creasing velocity ratio (uj/ Ug) for systems with approximately the same 
density ratio (p;/Po). The centerline decay is indicative of the rate at 
which the streams mix. If the centerline decay decreases, the rate of 
mixing decreases; also the length of the potential core increases as the 


*The region less than two nozzle diameters from the nozzle exit 
plane. 
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centerline decay decreases. The length of the potential core was deter- 
mined by the method suggested by Zakkay and others (Ref. 6). The 
centerline composition was plotted versus axial distance on logarithmic 
paper. A curve through the data was extrapolated until it intersected 
the line corresponding to 100-percent jet gas concentration. It was 
assumed that the intersection defines the core length. A similar method 
was used to determine the velocity core lengths. The velocity core 
lengths were longer than the corresponding composition core lengths, 
which indicates that the transport of mass is more rapid than the trans- 
port of momentum. 


It is interesting to note that the result obtained — that the mixing 
decreases as the velocity ratio increases — is exactly opposite to that 
predicted by the Prandtl eddy viscosity model presented in Section I. 


The density ratio of the inner jet gas to the outer stream gas was 
increased in Series II by heating the outer stream gas. The density 
ratio was increased even more in Series III by using air as the jet gas 
as well as the outer stream gas. Figure 7 presents a comparison of the 
centerline velocity decay for the three different density ratios at approxi- 
mately constant velocity ratio. The curves show that the mixing rate 
decreased with increasing density ratio for constant velocity ratio con- 
ditions. Since the mixing rate decreases with increasing velocity ratio 
and with increasing density ratio, it might appear that there would be a 
correlation between the mixing rate and the mass flux ratio (pjuj/ Poo). 


Another possible reason for considering this correlation is that the eddy 
viscosity model of Ferri, presented in Section I, is based on mass flux 
difference. This correlation was attempted using the centerline decay 
data, but it did not prove to be a valid correlation. In other words, the 
mixing rate did not decrease monotonically with increasing mass flux 
ratio. 


Zakkay and others (Ref. 6) found the composition decay downstream 
of the potential core to follow the relationship Ce = (x/xX9¢)7?. Figure 8 
shows a comparison between Zakkay's results and the data from this in- 
vestigation. By considering the expression C, = (x/x9c¢)"", where 
Zakkay found n = 2, the results from this investigation gave n = 1.7. 
Zakkay also presented a generalized expression for potential core length 
variation with mass flux ratio. Figure 9 shows a comparison between 
Zakkay's generalized expression and data from this investigation. The 
dashed curve represents the function x9/Ro = i3Va , and it agrees 
reasonably well with the experimental data. The slopes of the dashed 


curve and the solid curve from Zakkay's expression are the same, but 
the constants differ by 40 percent. Zakkay's generalized expressions 
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do not seem to apply in this higher velocity ratio regime; however, dif- 
ferences in the initial conditions, such as boundary-layer thickness, 
may account for the discrepancies, 


Figure 10 is a comparison of the centerline decay of composition, 
velocity, and total enthalpy for representative tests in Series I and Il. 
These curves show that mass and energy are transported more rapidly 
than momentum. 


The composition and total enthalpy decay is approximately equal 
for all of the tests conducted. This result may be expected in the 
Series I tests because the enthalpy was a stronger function of composi- 
tion than of temperature, and the temperature difference between the 
streams was only 100°R. But, in the Series I] tests, the temperature 
difference was 500°R, and the effect of different turbulent transport 
properties for heat and mass should be evident. Since nondimension- 
alized composition and total enthalpy were equal for the Series II tests, 
the turbulent transport coefficients appear to be equal. Consequently, 
the turbulent Lewis number is unity. 


3.4 PROFILE SHAPES 


The nondimensional composition for all of the hydrogen-air tests 
is plotted versus radial distance in Fig. 11. The composition is non- 
dimensionalized and normalized by dividing by the centerline composi- 
tion. The nondimensionalized radius is obtained by dividing by the 
radius at which the composition is one-half of its centerline value. A 
band representative of the data in Fig. 11 is presented in Fig. 12. 
Curves obtained from three different mathematical functions (cosine, 
exponential, and power law) are compared with the experimental data 
band. Each of the curves presented gives reasonably good agreement 
with the experimental data. 


Nondimensional velocity is plotted versus radial distance in Fig. 13. 
The parameters were nondimensionalized and normalized in a manner 
similar to those in Fig. 11. The solid curve, which fits the data well, 
is the cosine function. The data in Figs. 11 and 13 are representative 
of all hydrogen-air test conditions and all axial locations in the second 
regime. The curves indicate that profile similarity is an adequate 
assumption for engineering calculations in the velocity ratio range con- 
sidered. The data outside R/Rmc and R/Rmy = 2.5 are not considered 
because the data in this region have a small hydrogen concentration, a 
small velocity difference, and consequently, a large uncertainty in the 
calculated values (see Appendix III). The nondimensionalizing radii for 
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‘composition and velocity are presented in Table IT so that the actual 
radii may be obtained from the curves in Figs. 11 and 13. 


Figure 14 presents a comparison between typical experimental com- 
position and velocity profiles from the same axial location for two dif- 
ferent tests. The nondimensionalized composition and velocity are 
plotted versus radial distance. The composition profile is slightly 
wider than the velocity profile in each case. 


In the first regime, the statement of profile similarity is 
u-U, _ rR - R; 
Uy—- do ; ( b ) 
The term (Rj) is the radius of the potential core, and b is the mixing 
zone width. Since it is very difficult to determine Rj or b accurately, 


the expression 
R <. R.s 


which is related to 


is used in Fig. 15 to illustrate velocity profile similarity. The solid 
curve is obtained from the cosine expression, 


wy = th + cos a (Rp + 0.4)| 


and agrees quite well:with the data. 


3.5 MOMENTUM, ENERGY, AND MASS TRANSPORT 


In Fig. 16, the nondimensionalized velocity is plotted versus the 
nondimensionalized total enthalpy for all of the hydrogen-air tests. The 
velocity was nondimensionalized by using the jet velocity, instead of the 
centerline velocity, in this case. The dashed curve is the result for 
unity Prandtl number, and deviation from it is indicative of nonunity 
Prandtl number. These same parameters are plotted in Fig. 17 for the 
air-air test, and a distinctly different trend is evident. 


There appears to be a definite relationship between the transport of 
energy and the transport of momentum. However, it is not obvious how 
to: estimate the Prandtl number for this relationship. It is possible that 
a relationship between the velocity and total enthalpy, such as the ones 
presented in Figs. 16 and 17, may be more useful for engineering 
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calculations than the Prandtl number. This is even more probable if 
the Prandtl number must be treated as a variable. 


Composition is plotted versus nondimensionalized total enthalpy in 
Fig. 18. The solid curve indicates the result for unity Lewis number. 
The experimental points deviate only slightly from the curve, which in- 
dicates that unity Lewis number is a good assumption for these data. 


In all hydrogen-air tests, the turbulent transport of mass and total 
enthalpy is more rapid than the transport of momentum. This result is 
in agreement with other investigations reported (Refs. 3 and 4). 


The data presented in this report are being further reduced to deter- 
mine the turbulent transport coefficients. This is being done by numer- 
ically solving a set of integral equations involving density, velocity, 
pressure, composition, and enthalpy. The ratio of the transport coeffi- 
cients will be used to determine the variation in Prandtl number indicated 
in Figs. 16 and 17. The results of this work will be presented in a future 
report. 


3.6 BOUNDARY-LAYER THICKNESS 


Total pressure distributions downstream of the nozzle lip were 
measured to determine the boundary-layer thickness. The measure- 
ments were made with a stainless steel probe of 0.004-in. OD and 
0.001-in. wall. The total pressure traverses were made at axial sta- 
tions less than 0.010 in. downstream of the nozzle lip. Data for one 
hydrogen-air condition and two air-air conditions are presented in 
Fig. 19. The combined width of the inner and outer boundary layers (A) 
is approximately 14 percent of the inner nozzle radius. The thickness 
(A) is measured between the 99-percent velocity points in the inner and 
outer boundary layers and includes the thickness of the nozzle Lip. 


3.7 STATIC PRESSURE 
Typical static pressure distributions are presented in Figs. 20 
and 21. The pressure gradients are relatively large in the near field 


region; however, ina short distance, the gradients are smoothed to 
less than one-percent variation. 
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SECTION I¥ 
CONCLUSIONS 


The primary differences between this investigation and other in- 
vestigations of turbulent mixing are that (1) the initial boundary-layer 
effect was reduced by designing the nozzles to minimize the boundary- 
layer buildup at the entrance to the test section, and (2) the velocity 
ratio range was higher than that previously reported. The following 
conclusions may be drawn for subsonic, axisymmetric mixing, from 
this experimental investigation. 


“he 


For hydrogen-air mixing systems, the centerline decay of 
composition and velocity decreased with increasing velocity 
ratio (jet velocity/outer stream velocity). Also, the center- 
line decay of velocity decreased with increasing density ratio 


_ for systems with the same initial velocity ratio. 


Similarity of composition, velocity, and total enthalpy profiles 
is a valid assumption for engineering calculations in the velocity 
ratio range from 2.4 to 6.3 for hydrogen-air mixing systems. 
The commonly used expressions for the profile shapes, such as 
the cosine function, the three-halves power law, and the error 
curve, are representative of the shapes. 


There appears to be a definite and consistent relationship be- 
tween the transport of momentum and energy for hydrogen-air 
and air-air mixing systems. The Prandtl number is not unity, 
and it is not obvious that it is a constant. 


Unity Lewis number is a valid assumption in the velocity ratio ' 
range considered, at least for streams which have moderate 
temperature differences (jet gas temperature/outer stream tem- 
perature less than 1.0 but greater than 0.5). 


The data presented in this report are being further reduced to deter- 
mine the turbulent transport coefficients. The results of this work will 
be presented in a future report. 
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Fig. 7 Centerline Velocity Decay for Systems with the Same Initial Velocity Ratio 
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XIXoe 


Fig. 8 Centerline Composition Decay, Comparison with Zakkay’s Results 
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Fig. 11 Composite Second Regime Radial Composition Profile for All Hydrogen-Air Tests 
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Fig. 15 Composite First Regime Radial Velocity Profile for All Hydrogen-Air Tests 
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Fig. 16 Velocity-Enthalpy Relationship for the Hydrogen-Air Tests 
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Fig. 17 Velocity-Enthalpy Relotionship for the Air-Air Tests 
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Test 
Number 


IA 


IC 


. 2224 


2907 
3556 
4292 
5172 
6081 


1896 
2710 
3577 
4403 
5212 
6047 


1936 
2735 
3252 
3619 
4060 
4452 


1208 
1659 
2343 
2859 
3521 
4168 
4445 
4847 
5489 


. 2028 
. 2283 
. 2876 


3213 


. 3689 
. 4478 
.5337 


TABLE I! 
NONDIMENSIONALIZING RADII 


pee | 
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TABLE I] (Concluded) 


Test 


0.1911 0.0198 0. 0216 
0. 2460 0. 0202 0. 0203 
0, 2644 0.0210 0.0211 
0, 3473 0.0237 0. 0226 
0.4332 0. 0286 0.0264 
0.5158 0. 0320 0.0297 
0.5992 0. 0340 0.0318 
0.1705 0.0169 0.0196 
0.1890 0.0182 0.0196 
0.2161 0.0192 0.0192 
0. 2451 0.0215 0.0200 
0. 2735 0. 0231 0. 0208 
0.3128 0.0241 0.0217 
0. 3694 0. 0273 0.0243 
0.4191 0.0297 0.0258 
0.5178 0.0332 0.0282 
0. 2019 0.0172 0.0178 
0.2397 0.0195 0.0180 
0. 2906 0.0221 0.0187 
0.4377 0.0270 0.0232 
0.5178 0.0275 0.0228 
0. 0. 

0. 0. 

0. 0. 

0. 0. 

0. 0. 

0. 0. 

0. 0. 

0. 0. 

0. 0. 
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APPENDIX Ill 
RELATIVE ERROR ANALYSIS 


The measured parameters are listed in Table IJI-I with an estimate 
of the measuring accuracy. 


In the following error equations*, y is assumed to be a constant. 
This assumption is justified by considering that the maximum deviation 
of y is from 1.4 to 1.7 percent. This deviation was obtained when the 


outer airstream was heated to 1050°R. 


The relative error equation for ae is 


du | Cc aH (ip = Ho ti 0) itty 
. [e+a-o | Hh fs - a - ©) a) ¢ arr 2 Nis 0 Ha 
h 


The relative error equation for density is 


Mp Soh. Ye plik ee dP, (@o-@j) C_ dC 
p Ty y P, y Ps (@9 - W)C + w C 


The relative error equation for velocity is 


dus 1s dT, Tt; 5 a nae “4 (@9-@;)C dC y-1 1 T,] dP, 
—_— Ss = = + - V7——_———_ eee i - —_—_—_—_—_— — 
2 T, Tsy’ 2y (ow, -@j)C +@; C 2y yM? T,] P, 


u 


These equations are written for the mixing flow field conditions. 
They are applicable for the plenum conditions if the appropriate meas- 
ured quantities are used. 


The terms (Rmc, Rmu, and Rp) are selected from curve fits and 
do not depend on the accuracy of the measurement of a single point. It 
is estimated that they are accurate within +3 percent. 


The relative error equation for y, is 


dy ud ue due uy (u — ue) dug 
Yu u-U, u Ug ~ Um Uc {u — ug) (u, — ug) Ug 


*To obtain these equations, the equation for the parameter in ques- 
tion is differentiated, and the result is divided by the original expression. 
The equations for enthalpy, density, and velocity are contac’. in 


Appendix IV. 
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The parameter (y,) is susceptible to large inaccuracies when u ap- 
proaches ug near the outer edge of the mixing zone. For this reason, 
the uncertainty of the data for R/Rmy > +2.5 is very large, and the 

data are omitted from the curves presented. The same situation is 
true for the composition data. The low composition measurements, 
which have large measurement uncertainty, occur at the outer edge of 
the mixing zone. Consequently, the data for R/Rmc >+2.5 are omitted. 


TABLE Ill-i 
MEASUREMENT ACCURACY 


Estimate of 
Measurement Accuracy, 
Percent of 
Absolute Value 


Parameter 


All Considered 
All Considered 
All Considered 
All Considered 
All Considered 
All Considered 
All Considered 
All Considered 
All Considered 
>0.05 to 1.0 
>0. 02 to 0.05 
+50 >0.01 to 0.02 
Uncertain 0.0 to 0.01 
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APPENDIX IV 
CALCULATION PROCEDURE 


Gas mixture properties were calculated from the measured param- 
eters by using the relationships listed below. Tables of specific heat 
and enthalpy for air and hydrogen over a temperature range from 500 
to 1500°R were put into the computer program. The information was 
taken from gas tables (Ref. 10). 


Jet Mass Flow 


where K;, a dimensional proportionality constant, is furnished from the 
orifice calibration. 


Total Enthalpy H = CH (1 - C)H 
= h + = a 


where Hy and Ha are obtained from the enthalpy table for the measured 
total temperature. 


Q) Wo 


Cla@o - @j) + o, 
Specific Heat 
C, = C(Cpj) + (1 cn C) Cpo 


where Cp; and Coo are selected from the table of specific heats for the 
measured total temperature. 


Ratio of Specific Heats 


Mach Number 


Static Temperature 


Density 


Velocity 


u= N y(Ps/p) M? 
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APPENDIX V 
CONSISTENCY CHECK 


The hydrogen mass flow was calculated at each axial station where 
data were obtained. The mass flux was numerically integrated to a 
specified control volume. The expression for the calculation is 


R* 

Wo= 2n f puCRdR 
Numerous control volume sizes were tested. As long as the control 
volume was large enough to extend to a location where the hydrogen con- 
centration was negligible, its size had little effect on the value obtained 
from the integration. The ratio of the values obtained from the integra- 
tion (W) to the hydrogen flow measured by metering orifice (Wj) is shown 
in Figs. V-l and V-2. The results are reasonably good, considering 
that the hydrogen mass flow is only from 0.4 to 1.7 percent of the total 
mass flow in the mixing system. 


A possible reason for the low calculated mass flow is that the probe 
was not on the vertical centerline when the data were recorded. There 
was no fixed reference point from which the vertical position measure- 
ments were made; therefore, only vertical position data recorded during 
the same test run may be compared. The data for tests I A andI D were 
recorded during the same test run; the vertical position variations of 
tests J A andI D were 0.050 and 0.070 in., respectively. TestI A gave 
good agreement on the consistency check, but test I D gave low values of 
mass flow calculated by the integration technique. 


The following calculation was made to get an estimate of the effect 
of probe position error: The data from test 1 A were assumed to be on 
the true vertical centerline, and the integration was made for a traverse 
0.060 in. below the centerline. The calculated change attributable to the 
assumed probe alignment is from 2 to 5 percent, which is too small to 
account for the low values obtained during test I D. 


The tests in Series II showed a trend toward high values of inte- 
grated mass flow near the nozzle exit plane; then the values decreased 
with increasing axial distance. The reason for this result is not known. 


Even though the consistency check does give poor results for some 
tests on an absolute basis, the repeatability of the data is very good. All 
other data from tests which gave poor consistency checks agree well with 
the results for tests which gave good consistency checks. This fact is 
illustrated by the figures which present data for more than one test. The 
data appear to be better than the consistency check indicates. 
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APPENDIX VI 
EXPERIMENTAL DATA 


Nomenclature for Tabulated Experimental Data 


Cc Mass fraction of jet gas 


CP Specific heat at constant pressure 

CPJ Constant pressure specific heat of jet gas 
CPO Constant pressure specific heat of outer stream gas 
GAM Ratio of specific heats 

GAMJI Ratio of specific heats of jet gas 

GAMO Ratio of specific heats of outer stream gas 
H Total enthalpy 

HJ Total enthalpy of jet gas 

HO Total enthalpy of outer stream gas 

MW Molecular weight 

MWJ Molecular weight of jet gas 

PS Static pressure 

PT Total pressure 

PTJ Total pressure of jet gas 

PTO Total pressure of outer stream gas 

R Radial distance 

RHO Density 


RHOJ Density of jet gas 
RHOO Density of outer stream gas 


RS Control volume radius 

TS Static temperature 

TSJ Static temperature of jet gas 

TSO Static temperature of outer stream gas 
TT Total temperature 

TTI Total temperature of jet gas 

TTO Total temperature of outer stream gas 
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Velocity 

Velocity of jet gas 

Velocity of outer stream gas 

Jet mass flow as measured by a choked orifice 


Axial distance 
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TEST NUMBER -1C 
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TEST NUMBER - 1 E 
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TEST NUMBER -I1A 


Jet Gas — Hydrogen 
Outer Stream Gas — Air 
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20 7038-01 0 2.720F-02 
W.9dbE-01 2,04RE OL 2.0948 09 DAIL -O2 
2.S10F-O1 2.0348 Of 2.0998 03 VoMabC-02 
2.919808 2.676€ OF 2.0918 09 E,R2TF-02 
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2,5098-01 2ZoATOE OF 2.0008 OF €,823€-02 
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asietuset Tediet 190101 mole?usen TroL en oTOLgss) cpoletuse @1 = Gann OndopeasF? Sh 
Oo OOL RE O2 E,O3878 GS F,.900RT G2 O.d249E O2 1.0002 OF Z.O2F1E OF 2.490801 1.3778E OO 3,3910P-02 
olzer € vse Bi an Fed F2etFl GHO(OMsET BQ ULPFFEECE mMIidtUseL 
Oo 130E-O2 3.3KE-04 @OLSE-01 1. 370E OO 2.3968 OF 2,050€ 0 O2 3.O22E-O2 T7898 O2 2.9292 OF 
VoMEE-O2 BW LIZE-04 2, 910R-O1 F.2PSE OF 2,0848 OL 2.0590 @2 BeAATEOP P.PAPF OF 2.9912 OF 
SoPlOE-O2 2,080E-04 2,500F-01 1.2708 00 Ot 2.0398 O2 3,9008-O2 T.TEZE OF 2.9228 
Vee DPOF 1.000803 2,E298-03 1od7OF 00 OL 266938 @2 3.030802 PoTSZE O2 2649258 
HoMIE-OF «Ww IIE-O) 2 3OPE-O8 8. OLE OO Ob 2.0308 OF E,O2ZE-OF F.O3sF OF 2,939¢ 
HeBQZEO? WeHOPEOZ PW SHCE-Ol 2.2008 00 Ob 2.0318 OF B.CPIEOZ 7.7668 OF 2, AGE 
6.037202 b,1FE-O2 2, 7OTE-O1 8.3948 00 Ob 2,040€ OF F.1VOE-OF SLATE OF 22,0307 
CLATE-O2 9%, bame-O2 4, 300F-O1 1297 2.0038 O2 2,038€-02 3738 OF 3.9108 
VWSISE-OF 1. ISTE-OL 6, FOPF-O1 € 2, 000€ O2 2.022F-02 1.061F OF 4.4968 
VeTEVE-O2 1,RbOEoON *?0F-OL 2.037€ O2 L.OO8E-O2 1,2138 OF S.ASZE 
G.427E-03 2 PWE-O1 9, EFOI-01 2.O30€ O2 2.466€-02 1.260F OF 6.4828 
2.001504 2,4017-01 1.0118 00 2.0328 O2 1.O4bT-O2 1.9ROE OB 6.6198 
6.109802 B35 7TE-O) 6, VO8E-O8 2.0338 O2 UeASKE-O2 L.dTOE OF HL. DAGF 
UIETEO? 2.03901 B,OAIE-OL 1.9908 2.0%€ O02 Ue568E-O2 2.2308 OF 3,6308 
eR 1GTESO? L,S7LE-O 7, 499E-01 2.0388 O2 L.360F-GC2 LLLhOe OF 3.0078 
HP. POAT-O2 ULLOKE-O) DL OTLE-O1 bolPE 00 Ol 2.0026 OF 2.27102 FAME OF 4, 2482 
THOHTE“O2 O,B49E-O2 4,440F-08 1.999 00 OL 2.0078 O2 2. TIBF-OF &oO39F O2 S,0dKE 
HHAAIVEO2 2,098EO2 BA 26OF 01 163908 00 OL 2.050€ O2 3,395€-02 7,092 O02 2,030€ 
*S.229-02 9.977E-99 2.0TSF-O1 1.ddEF 00 Ob 2,033t O2 @WF-OF 1%, MBE O2 72,9678 
6.660802 2.91GE-Ob bLdTSE OO 2.8708 OL 2.0988 O02 3.83720-0F T.AOOE O02 72,9408 
-€, $0102 2.40201 1.3708 00 2 € ot 2.091€ O2 Y.MAGE-O2 7. 7SOF OF 2.5708 
0. 3306-02 2oSIZE-O1 0.dTAE 08 € ob 2.051 OF 3.449F-O2 T.799€ OF 2.52K0 
7 MbLal OR 2.49E-O1 1.376 00 © ot 2.09t€ OF O2 3.924F-O2 FeOLZE OF 2.4 7HE 
whe L22t-01 2.40G-RI 0.3028 00 2.002€ Of 2,031€ O2 4.063€-02 1.2638 O2 2.9018 
o1pes @evezs UsLOASEACL Musionsmmel § FEIIOS Peswprel §=6ceateeuse Ol = GAMD Remsiense? eb 
bo LPNSE-OL O,PDISEOR 1.67SOE-O9 2.01086 O08 0.01938 O2 2.33042 CF B.43198 00 1.40268 68 8. 98 NE-03 
matetusen Oosnen 220101 +O 1etuset Troves PTOLPSER «<CealeTuse BL 8 6CAnn aHaoiease2 $B 
1.09952 OF B.4005E OP 9.00002 O2 2.22542 OF 1. O4G4E OF 2.62602 OF 2.49008-01 1.37702 OO 3,04708-02 
Bieri CHietuse Ob = GAA ESLFERD E9¢O) «= OrmrtemsPy 31 ULEESSECD NierUseL 
DIR 02 Ob 2.0662 03 49.9308 OF 3,69F-02 7,597F O2 2.6238 OF 
TobL 2-02 Ob 2.0908 eerie: VePS9E O2 2.5068 OF 
1812-02 Ob F,O00E v.07 2.90398 32 
DAA EOF Ob €.0428 JoEO4E-O2 23198 OF 
Ob 2.0528 30883002 2.4068 02 
i) Ot €.OSt€ DOLZE-OZ TVoMSE OF G.9408 OF 
VOLO? Le zTIF-O2 ot £.0908 V.648E-O2 ToPOLE O2 2.6228 OF 
4.3290-02 E,3DA0-02 E.207E-C1 1.3907 OL €,0098 B.432CO2 7.907 O2 2,709" O2 
DeVHME-O2 5.290802 6,192F-08 Ob 2.0408 OF 7.608 O2 2.630F-02 d 322408 02 
PALOF-O2 Fe OO2E-02 4,6805-O1 Ob €,0438 OF V.265E OF 2,.670F-02 3.3908 02 
2.780€-02 1,0998-01 Ol €.0622 03 6.970 O2 2.279€-02 Se 179 OF 
2.2308-02 L.E1ZE-08 Ob 2.000F OF e. Tet O2 Set-o2 4.278 02 
V.O20€-02 be 42TE-00 6.3238 00 2.0388 03 €-02 4.009€ OF 
b.3*re-02 @.070E OO 2.0372 OF B.d7CE OF LeNLVE-OF 8.4672 OF 3.2272 OF 
: 28-03 0.9798 00 2.08 03 U.7O0E-02 1.2108 OF 
478-0) @. 3062 00 2.0368 03 1.2308 OF 
3. 4008-0) 8 00 2.039 OF €.4938 O2 b. 260 OF 
9.19404 1.900F-0) 0.300F-08 1.3008 00 E.O33E OF &.004F OF 
*2.3308-03 L.9O2E<OL O,S13E-01 1.3098 00 2.0388 03 6.0986 OF 1. 68EE-02 
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7.24007 1.000F-09 00 2.0508 08 3.867E-02 E.SI4E OF 
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TeQUIEOL b.O00F-23 2,40 7F-Gt 1.3088 00 2.0698 OF 8,907 O2 4.28LE-O2 6.858E 02 


76 


USIP E/SOC) 
0.3760 09 


UBNET/S2CL 
Tob2O7F O2 


Tre) 
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Lal de) 
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AEDC.-TR-68-133 


TEST NUMBER - II C 


Jet Gas — Hydrogen 
Outer Stream Gas — Air 


opel Agie2> WILOMSSECL Anyi OmsNMEE =6—T22IAt orsieset  cesretuse 61  GAng CeBIT ORF Ba wsie2/90C8 
ToL WEG? 44h DWH-OR LoGITI-Oy E.OLGOF OO 6.10008 OF 2Z.DIITE OD F.AdETE OO §.401SE OH 4.63368-03 1. TIOCE OD 
wslatusel Tastal 1$0¢al c@iarusad rrolal OTe Pare Crareftuse oF t canted RHOOL OMsFZ DL VOLPE /S ECE 


2. OF63 OF 6, PLOOE OZ §.O1608 OF 2,59808 O2 L.OFOPT OD 2.4346F OD 2.SO3OH-OL LL ITIOE 08 3.86108-02 7.33008 OF 


anoleney Bt Ute V/BIC) wietusel “¥tieh —FYewy 
Od D.M40E-02 1.5568 O2 2.90TE O2 1605R8 OB 2004618 OD 
ToGL2E OF BW SOLE OP §.06CH OD 2.4088 OD 
T.419€ OF 2.53968 OF LeOSHE OD 2.4428 OF 
ToLTIE OZ BZ446E O2 LOL4E OF 20408 OD 
_1, 7308. 03 oF O2 2.3348 _ 0) 


aiete (4 colartuse at cu La) 
Q.OSLE“O2 2.120E+0O4 2,S10F-O1 ILD7TE OO J2.MAVE O1 
O.FIE-O2 1,0002-05 2,506F-O8 LedITE OO 2.8078 OF 
$.S10F 02 1.0008-05 2.5046-01 1.3772 OO 2 ot 
D.O92E-O2 L.O00E-0S 2.490F-91 §.380F CO 2 ot 
he SO TE*O2 _ MoE NTL-O)__ Do IAS 09. 1.4028 OO 2604 76 09 26 . 
WE-DL B41OE OO 1.407E 2.0286 OO 2.0942 03 


226-03 


Le? 28-02 oz ULE-Od) LoVITE OD 03 0108 02 2.9908 6 
OAS 1-03 BOL 3.4282 00 1.4028 be OO 2.09E 03 oz 4%E-0O3 $7828 O38 03 OdE C2 3.dd4E OF 
“OL 3.4%E OO 1.8018 00 OO 2,.093F O8 oz 2BOTE-0} 1. 7OSE OD 0) $2 02 2.3348 03 


oo OME OF 3.dd4E OD 
OD G.2ONE O2 BZoddIE OD 
O2 7488 O2 2.4278 33 
O02 1.058 Od 2.4626" 

OF L.O8TE-O3 2,445€ OD 
O2 b.000€ OF 2.4558 03 


448-01 3,228F 00 b.402E 00 
Y2E-01 2.8108 06 1.603€ 00 
o275E*O9_ 20 5206-0) 8 .292€ 00. 
2.2ISF-O5 =—-2,FO4E-O1 = LWDTTE 00 
L,0008-03 2,9077-O1 1.37SE 00 
1.Q00E-05 2.00¢F-01 1.377E 00 


OO 2.0976 OF S.913E O2 4.°276-03 1.7206 03 
00 2.09F O} G.O1dE O2 3.567E-03 1.6262 OD 
OL_2,09*E OB _@. 3006 O2 4,1236-02 0,999E 02 
Ol 2.0°)€ O} L.O1SE O3 3.800f-02 7.0088 02 
OL 2.009 OD} 1.022E OD 3.8502-02 7.5198 OF 
Ol 2.0R1F CS 1.019 OF 9,830F-02 T.A9SF 02 
beOWVEOS 2.5048-01 LedTTE 00 O1 2.0772 03 1.08SF OD 1.8952 02 02 10006 Od 3.4418 OD 

beQO0E-04 2.502t-01 L.377E 08 2.097E Of 2.O7TE O} 1,007E OF 3.0686-02 T.SASE O02 O02 L.OSDE Od 266438 OD 
—The262k-O1 LeOOOF-O5 2SOMF-O1 _1oITAE OO _2.097E OL 2.0771 OD} 1.006F OD D.RTRE-OF T.S92E O2 2.5958 O2 1.0S0E OD 2.4438 OD 


oees OTs VITON/TECT MejlOn/MOLEL © tsiel PEILOSFE §«CHIIOTUsO BL 8=6—Gamy CML OM /PY BE UTE 2 tect 
. is rete SB, TO-Oy 1.62 T8E-O2 “reteee 02 6.000RE G2 2.32018 OD D.ADOIE OO 1.00868 OO 4.90008-03 1.91008 09 
wMgiatused = «Fb aNe Vota wolMTusal veqial PEOLPSFL §=«CPRI@TUsN BE «66ORO )8=— BHMOLOR/PY DD UOLOE SEC! 
2OS4SE OF 9.0070E O2 1,00%F OF B.DSODE O2 8. OSO3E OF 2Z.4DHE OD 3.502%E-O1 1.37728 CO D.ODOSE~O2 T.84O2f OD 
auert c Ceiefuse wt (OAM ne OL t Vsiat meoiewVey th gS py ae Wiel vrvwr 
0.090E-02 140-03 D.SSTE-OL 1.37TTE OO 2.0618 © 2.0008 03} 1.0108 ©3 3.7S0E-O2 1.0408 2.005€ O2 1.0308 OD 2.441E OD 


SelB GEO? be2IHOF—OD 2.54601 LedDTTE CO 2.080F OL 2.0078 OF 1,0098 O3 J.7OHE-02 7.0208 +i 2.395€ RA LeOSTE OF 2.048 OD 
SoM62E-O2 = LALGE“Od = 2.9388-O1 1.3781 00 2.0626 O} 1.0068 OF B.7TATE-G2 7.82008 OZ B.5AIE OP 1,054E OD BVOOE OD 
GALTIE“O2 2.0528-03 2,567F-O1 LodTOE 00 2.0838 OF L.OCOLE OD F.TOUL-O2 %,809F O2 2.5°9E O02 
— oe 24bl 02 GaP ITEOD 2 PHIE-O1_ _badSHE_ 00 8.0392 O2 3,.0256-02 F.370E G2 2.532E O2 
W.60TE-O2 LeSOd1-O) 7F.1GeF-O1 1.,300F 00 O.H16E O2 L.978E-O2 @,0018 O2 4.7808 OF 
WOM O® AeLBAE-O1 8.5028 OO 1.4018 00 OoO7E O2 F.PA4E-O} 1.352E OF D.GTIE O2 
VedeeloO2 8.334801 2.004F OO 1.402F 00 SeASLE-O} L.BOLE OF 8.T22F OD 


00 2.058 
00 2.0506 03 


&.905F-O} 9,.509F-0} B.292r OO 1.402F 00 00 2.050€ O38 e.7S0E-O3 1 1.9001 03 
#2,409F-03 8 MEF-O1 J.IMTE CO 1.4018 00 OO 2.053€ 6» 4.5%1-03 1 2.0018 03 
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